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ABSTRACT 

DAOSPEC is a Fortran code for measuring equivalent widths of absorption lines in stellar 
spectra with minimal human involvement. It works with standard FITS format files and it is 
designed for use with high resolution (R>15000) and high signal-to-noise-ratio (S/N>30) spectra 
that have been binned on a linear wavelength scale. First, we review the analysis procedures that 
are usually employed in the literature. Next, we discuss the principles underlying DAOSPEC 
and point out similarities and differences with respect to conventional measurement techniques. 
Then experiments with artificial and real spectra are discussed to illustrate the capabilities and 
limitations of DAOSPEC, with special attention given to the issues of continuum placement; ra- 
dial velocities; and the effects of strong lines and line crowding. Finally, quantitative comparisons 
with other codes and with results from the literature are also presented. 

Subject headings: Data Analysis and Techniques 



1. Introduction 

Most of the information we acquire about the 
chemical compositions and atmospheric parame- 
ters of stars comes from spectral absorption fea- 
tures, both molecular bands and atomic lines. 
One of the most powerful and well-understood 
techniques to disentangle the chemical abundance 
from, for instance, the effects of temperature and 
gravity on absorption features employs measure- 
ments of equivalent widths (EWs). 



1 Based on data obtained from UVESOVLT, GI- 
RAFFEOVLT, HARPSOVLT located in the ESO obser- 
vatory of Cerro Paranal, Cile. Also based on data from 
FOCESOCAHA 2.2m telescope in Calar Alto, Spain. Data 
from the VALD and GEISA online databases has been em- 
ployed as well. 



Usually, chemical abundances are best de- 
rived when the EW analysis is based upon high- 
resolution (R= A/5A> 15000) spectra with rel- 
atively high values of the signal-to-noise ratio 
(S/N>30). This level of data quality was typ- 
ical of the best that could be achieved in the 
1970s and now it has become a sort of lower limit 
for the standard 2 quality expected for deriving 
accurate chemical abundances from optical spec- 
tra (A ~ 3200 - 10000 A). Recently, it has become 
possible to achieve similar data quality (R~10000, 
S/N~30) in the near infrared wavelength range 
(up to if-band wavelengths, A ~ 2.2/^m). 

As technology has advanced rapidly in the last 



2 Typically, R~40000 and S/N~50-100 are used in the opti- 
cal. 



few decades, telescopes of increasingly larger aper- 
ture have become available, as well as spectro- 
graphs that allow us not only to reach very good 
resolution and flux fidelity, but also to observe 
many objects simultaneously. It has now be- 
come commonplace to work on sets of hundreds of 
high-resolution spectra from fiber spectrographs 
such as, e.g., FLAMES mounted at the Euro- 
pean Southern Observatory Very Large Telescope 
in Chile (Pasquini et al. 2002). Our knowledge of 
atomic constants for many of the absorption lines 
present in stellar atmospheres has also advanced 
and we are now typically able to use several tens of 
lines for a few dozens of different chemical species. 

The actual process of measuring EWs has un- 
til recent years mostly been tackled with the help 
of interactive graphical software, such as IRAF 3 
or MIDAS (Warmels 1992). However, this proce- 
dure has two main drawbacks: (i) lines are usually 
identified by the user in a non-automatic way, and 
the measurement of each individual line requires — 
to say the least — several seconds of direct human 
involvement (lack of performance); (ii) a certain 
degree of subjectivity is involved in the measure- 
ment process, particularly due to the continuum 
placement by the user around each line; this means 
that different users could produce different results 
from the same data (lack of robustness). It may 
also be noted that the tedium of the process itself 
often discourages a thorough testing of the proce- 
dures via repeated measuring under different sets 
of assumptions and protocols. 

Therefore, given the vastly increased volume 
of data and the better data quality available 
these days, the community has started writing 
(semi-)automatic codes that improve both the 
performance of the measuring process and the 
robustness of the results. Most of these codes are 
not really public, in the sense that they are de- 
signed to meet the specific goals of some research 
group and are provided on a personal basis to 
close collaborators. 

DAOSPEC is a Fortan code to measure EWs 
of atmospheric absorption lines in high resolution 
spectra. Our goal is to provide the performance 



3 Image Reduction and Analysis Facility. IRAF is dis- 
tributed by the National Optical Astronomy Observato- 
ries, which is operated by the association of Universities 
for Rcasearch in Astronomy, Inc., under contract with the 
National Science Foundation 



required by the increasing amount of data avail- 
able; the convenience and repeatability of a fully 
automatic routine; and availability to the whole 
astronomical community. With the help of sim- 
ple shell scripts, detailed in the DAOSPEC Cook- 
book (Pancino & Stetson 2008), a large number of 
stars with different characteristics and configura- 
tions can be processed in a fully automatic and un- 
supervised way. Effort has been put into making 
the code compatible with standard astronomical 
image formats such as the FITS standard (Wells 
et al. 1981; Calabretta & Greisen 2002; Greisen & 
Calabretta 2002; Greisen et al. 2006) and to use 
standard Fortran coding and libraries. 

One advantage of a fully automatic code that 
is also publicly available is that the results are re- 
producible and testable in variety of conditions of 
use and by many different users. Indeed, such a 
set of independent tests is presently under way in 
the small but growing community of DAOSPEC 
users (see also Section 4). 

In Section 2 below we describe the EW mea- 
surement process in general, and its implementa- 
tion in a few public software packages; in Section 3 
we give a general description of DAOSPEC itself; 
in the various subsections, we discuss in depth a 
few topics which are relevant to the most com- 
mon scientific cases, including radial-velocity de- 
termination (Section 3.1); continuum placement 
(Section 3.2); strong lines and blended lines (Sec- 
tion 3.3); uncertainty estimates (Section 3.4); and 
dependence of the results on the choice of input 
parameters (Section 3.5) and on spectral quality 
(Section 3.6). Some performance considerations 
are also presented (Section 3.7). Comparisons 
with other codes and with previously published 
measurements are presented and discussed in Sec- 
tion 4. Finally, in Section 5 we summarize our 
conclusions. 

2. EW Measurements in General 

Figure 1 illustrates how the EW of a spectral 
line is defined, starting from the continuum level 
(Ic ) and the specific flux at the line center (I ). 
The EW represents the strength of an absorption 
line, and it is defined as the width, in wavelength 
units, that a rectangular stripe of height Ic , in in- 
tensity units, would have if it had the same area 
as the actual line or — in other words — the wave- 
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Ic t [1 - fl(Ai)] 



(3) 








Wavelength (A) 

Fig. 1 . — Classical definition of EW: the width of 
a rectangular line (right) with the same area as the 
actual absorption line (left). See text for details. 
Ic is the continuum level while I is the specific 
flux at the line center Ao, measured upward from 
zero. 



length interval that would be covered by a hypo- 
thetical, perfectly opaque absorption feature that 
removes the same amount of energy from the con- 
tinuum flux. More practically, for an isolated spec- 
tral line in a discretely-sampled spectrum from a 
real instrument: 



EW = A\J2 



I Ci ~ h 



(1) 



where AA is the (constant) pixel size, Ic t is the 
continuum level at the wavelength of the i-th pixel, 
and li is the actual flux received by the i-th pixel. 

In a real spectrum, however, noise, spectral de- 
fects and neighboring lines can affect the perceived 
profile of an absorption line, and an EW is not nec- 
essarily readily measured ("IS tl direct summation 
of the pixel fluxes. Since the line profile is gen- 
erally dominated by the point-spread function of 
the slit image, which can often be approximated 
by a Gaussian profile 4 (see Section 3.3), often a 
numerical fit of a Gaussian-shaped function g(X) 
is adopted instead of the direct numerical summa- 
tion: 



g{X)=Ae-^- x ^l 2a \A>0 



(2) 



4 Different line profiles (i.e., Voigt, Lorentz) can be adopted 
for lines that significantly deviate from the Gaussian form. 



(i.e., the specific flux 7j assumes a minimum value 
Io = Ic, — A at Aj = Ao, and Jj — > Ic t for Xi far 
from Ao). The estimated EW thus becomes 



EW = Y J 



Ic, 



-AA 



g(X)X (4) 



where AA is the step of the integration process 
that can be made very small (— » dX). 

No measurement means much without an asso- 
ciated uncertainty. When fitting lines to obtain 
EWs with Equation 4 in a digital computer, a 
standard error can be computed for each line with 
various mathematical techniques that estimate the 
goodness of the fit (e.g., least-squares formal er- 
rors). These confidence intervals can also take 
into account the quality of the spectrum (R, S/N, 
and — to a certain extent — crowding and spectral 
defects, etc.). Cayrel (1988) gives approximate 
formulae for estimating the uncertainty of an EW 
as a function of spectral quality: 



._„. , V AA • FWHM 1 VAA • EW 
SEW ~ 1.5 — ~ 1.6- 



S/N 



S/N 



(5) 

where FWHM is the the full width at half max- 
imum of the line. (For a perfect Gaussian pro- 
file, FWHM=2.355 a.) Of course, as noted in the 
paper by Cayrel (1988), this is a lower limit to 
the actual uncertainty because the true continuum 
level is never perfectly known, but it too has an 
associated uncertainty that should be propagated 
through the whole calculation. However, the ef- 
fect of continuum placement (see below) is very 
difficult to take completely into account. 

Ways to estimate the uncertainty due to the 
continuum placement in each pixel have been ob- 
tained by Ramirez et al. (2001), Bohlin et al. 
(1983) and Levshakov et al. (2002) among others. 
The idea is that the derivative of Equation 4 has to 
be computed not only with respect to li but also 
with respect to the continuum level in each pixel, 
Id , and the total uncertainty then becomes 5 



5 For those programs that evaluate the continuum locally in 
two windows beside each line, the full uncertainty can also 
be estimated with the alternate method provided by Bohlin 
et al. (1983). 
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where Sli and Sid are the uncertainties on 7j 
and I Ci , respectively; e.g. these would be the 
Poisson errors \JTi and y/Id if the spectra are 
in units of raw detected photoelectrons. Other- 
wise, of course, the errors must be appropriately 
adjusted for any renormalizing of the data. As 
pointed out by Levshakov et al. (2002), the above 
formula still provides only an estimate of SEW be- 
cause when the spectra are resamplcd to constant 
wavelength bins, the errors in the single pixels be- 
come correlated. The correlation term that should 
be added to Equation 6 appears to be small, and 
as we will see in Section 3.2, there are other uncer- 
tainty sources that dominate the continuum error 
budget, such as line crowding and blanketing by 
very weak lines. 

A very simple experiment shows the impact of 
continuum placement on the EW measurements. 
An observed spectrum of star 2129 in the open 
cluster Collinder 110 (Pancino ct al., in prepara- 
tion) was normalized to an arbitrary "best" con- 
tinuum in the range 5500-6500 A and the EW 
of ^50 lines were measured with the task splot in 
IRAF. This spectrum has a resolution of R~30000 
and a S/N ratio of approximately 75 (per pixel). 
Later, the lines' EWs were remeasured adopting a 
continuum level higher or lower than the adopted 
"best" continuum by 1%, 3%, 5% and 10%, re- 
spectively. The resulting AEW values are shown 
in the four panels of Figure 2. As can be seen, 
except for a few lines with the smallest EWs, all 
measured line strengths are altered by approxi- 
mately constant values: <~ ±3-3.5 mA in the case 
of a ±1% misplacement, ~ ±9.5 mA in the ±3% 
case, ~ ±16 mA in the ±5% case, and ~ ±30- 
32 mA in the pessimistic case of a 10% misplace- 
ment. Of course, noise, resolution, line crowding 
and line blanketing also have an impact on the ac- 
tual AEW values, but the above experiment gives 
a quantitative idea of what happens when the con- 
tinuum level is badly misestimated (see also Sec- 
tion 3.4.3). 

In the following subsections we examine the var- 
ious operational steps that must be undertaken 
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Fig. 2. — Differences between EW measured on 
the adopted "best" continuum level and the EW 
of the same lines measured on a continuum level 
altered by ±1% (Panel a), ±3% (Panel b), ±5% 
(panel c) and ±10% (Panel d). The spectrum is 
that of star 2129 in the open cluster Cr 110 (Pan- 
cino et al., in preparation) 



when measuring the EW of a line. We will provide 
practical examples based on a few publicly avail- 
able programs. The IRAF task noao.oned. splot is 
taken as a good representation of those graphical 
interactive packages that allow the user to man- 
ually 6 measure EWs from their spectra; we be- 
lieve that inferences drawn from the IRAF pack- 
age will apply to other analogous routines, such 
as that provided within MIDAS, and others. The 
Fortran program SPECTRE 7 (Fitzpatrick & Sne- 
den 1987) is taken as an example of a more au- 
tomated routine and since it is part of the widely 



"When dealing with many spectra having similar charac- 
teristics, many users develop their own pipelines in the 
IRAF script language or the MIDAS equivalent, but these 
pipelines cannot really be considered here, since they are 
often highly specific to the problem treated and are seldom 
made public. 

7 http://vcrdi. as. utexas.edu/spectre. html 
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used MOOG package 8 (Sneden 1973) it serves as 
a very good reference in the following discussion. 
ARES (Sousa et al. 2007) and EWDET (Ramirez 
et al. 2001) are two programs written in CH — h and 
Fortran respectively. Although not entirely ob- 
jective (due to user involvement in the continuum 
placement, see Section 2.2), to our knowledge they 
are the only fully automated codes publicly avail- 
able 9 . 

2.1. Line Finding and Identification 

Within interactive tasks such as splot or SPEC- 
TRE, the processes of line finding and identifica- 
tion are generally joint and happen in the same 
initial phase of the measurement. By line find- 
ing we mean the process of specifying where in 
the observed spectrum are perceived dips repre- 
senting good candidate absorption lines. By line 
identification, on the other hand, we mean the pro- 
cess of specifying which of the perceived absorp- 
tion features corresponds to a particular atomic 
transition. The latter task requires some knowl- 
edge of the radial velocity of the star to infer the 
rest wavelength of the line for comparison with 
laboratory line lists. 

Within IRAF, splot requires the user to mark 
with the cursor, on a radial-velocity-corrected 
spectrum, the continuum position on both sides 
of a line. Apart from helping to define the con- 
tinuum level, this also tells the program the line's 
approximate location (initial finding). The pro- 
gram then fits the appropriate function and de- 
termines a precise wavelength for the line center 
(refined finding), as well as the EW. Finally the 
user confirms whether the measured line center 
is acceptably close to the laboratory wavelength 
(identification). 

SPECTRE also needs a radial- velocity-corrected 
spectrum as input, and accepts entry of the lab- 
oratory wavelength of a line (or reads a list of 
wavelengths from a file). The measuring process 
per se has a similar structure to that of IRAF: 
initial finding based upon the rest wavelength, 



http://verdi.as.utexas.edu/moog.html 

'We are aware of the existence of other codes that are not 
really public, and therefore are not considered in the follow- 
ing discussion. For instance, one such program, ROSA, has 
been written by R. Gratton (E. Carretta, private commu- 
nication) and another, fitline, by P. Frangois (P. Frangois, 
private communication). 



refined finding (and fitting) of the perceived ab- 
sorption feature closest to the rest wavelength and 
final identification by decision of the user. In the 
easiest cases, users just need to provide SPECTRE 
with the spectrum and a line list, and to check sub- 
sequently whether the program has done a good 
job. Thus, the procedure can be much faster than 
with IRAF. However, SPECTRE also provides 
several interactive tools to improve the fit in more 
difficult cases, such as continuum adjustments, the 
use of partial line profiles, and others. 

In ARES and EWDET, line finding is com- 
pletely automatic and independent from the iden- 
tification procedure. 

ARES uses local maxima of the second derivative — 
defined as zeros of the third derivative — to iden- 
tify the centers of absorption lines. In the paper 
that presents the code, Sousa et al. (2007) have 
shown (see, e.g., their Fig. 3; see also Section 2.3 
below) that this approach can enhance the detec- 
tion of highly blended lines not readily apparent 
in the original spectrum. However, the effect of 
noise dramatically propagates in the numerical 
derivatives even for spectra with high S/N ratios. 
Therefore, in practice the observed spectrum is 
artificially smoothed to reduce confusion between 
noise and real features before computing the nu- 
merical derivatives. This is especially important 
since the FWHM is left as a free parameter (see 
Section 2.3) and cannot therefore be used to dis- 
tinguish between an absorption line and a noise 
artifact. To avoid further problems, the input 
parameter lineresol tells the code the minimum 
distance between lines in Angstroms. ARES re- 
quires that the spectrum be pre-corrected for the 
radial velocity to identify detected lines with en- 
tries in the input line list on the basis of their rest 
wavelengths. 

In EWDET, the spectrum is examined, pixel by 
pixel, to find all the points that deviate from a pre- 
viously found continuum by more than 2a. Those 
points are further examined to see whether they 
have neighbors that also deviate from the contin- 
uum, and if so, the line finding is considered com- 
plete. The line identification is left to the user in 
the case of EWDET, and, to facilitate this task, 
the radial velocity is used as an input to compute 
the rest wavelength of each line. 

One side effect of all these methods, which can 
considerably slow down the procedure, is that the 
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radial velocity of each star must be known a priori. 
In some cases, as we have seen, the input spectrum 
must already be shifted to zero radial velocity. 

2.2. Continuum Placement 

Continuum placement is a source of uncertainty 
for any EW measurement, and there is as yet no 
generally accepted best practice in estimating the 
continuum level in a stellar spectrum densely pop- 
ulated with absorption lines. For instance, when 
writing a program for extracting EWs, one must 
consider whether to use a local or global contin- 
uum normalization, or a combination of the two. 
The approach employed in DAOSPEC will be ex- 
amined in depth in Sections 3 and 3.2. Here, we 
describe the procedures adopted by SPECTRE, 
splot, ARES and EWDET. 

SPECTRE requires spectra that are already 
roughly normalized to their continua. This is nec- 
essary because the user can only, if needed, re- 
adjust the local continuum level and not its global 
shape. Therefore it is optimal to normalize spectra 
that have steep or strongly curved continua before 
feeding them to SPECTRE. 

IRAF/ splot allows the user to mark separate 
left and right continuum levels relative to the cho- 
sen line, and these can have very different val- 
ues. However, since the continuum level within 
the absorption line is evaluated by interpolating a 
straight line between the two continuum samples, 
a preliminary continuum correction is again desir- 
able in the case of strongly curved or steep spectra. 
At least one of the private packages of which we 
are aware (E. Carretta, private communication) 
also chooses to normalize the continuum based on 
windows on the short- and long-wavelength sides 
of each absorption line. So, this package, like both 
SPECTRE 10 and splot, therefore requires (or gives 
better results after) a global continuum normaliza- 
tion, i.e., the procedure of dividing the whole spec- 
trum by a fitting function (usually a polynomial 
or a spline function), followed by minor localized 
continuum adjustments. 

ARES also finds the continuum level locally 
around each line, using a no~ rejection algorithm. 
The parameters of the fit and of the kg rejec- 



'We point out that SPECRE provides additional tools that 
help continuum fitting, such as the one colloquially known 
as boinking. 



tion, however, are highly customizable and flex- 
ible. This is the only point in ARES where the 
user has a strong impact on the outcome: trying 
many different combinations of parameters is im- 
portant for guaranteeing a good result, and an in- 
experienced user can obtain worse results than an 
experienced or extremely careful user, as can also 
happen with noao. oned. continuum. Moreover, the 
best combination of parameters will change for dif- 
ferent sets of spectra, and various users will have 
different opinions on how the continuum should be 
placed. In summary, continuum placement is the 
main reason why the results obtained with ARES 
cannot be considered entirely objective. 

EWDET chooses a global continuum normal- 
ization approach: the continuum determination is 
done by automatically fitting a curve to the whole 
spectrum with several iterations of point rejection 
above and below sigma levels specified by the user. 
This is not significantly different from what is done 
in the IRAF task noao. oned. continuum, and as in 
the case of ARES, it introduces a certain degree 
of subjectivity in the process. EWDET however, 
unlike other programs, stores in its memory the 
uncertainty of the continuum placement and uses 
it in the computation of the final EW uncertainty 
(see Section 2.4). 

The choice of the continuum normalization 
method is crucial and will be discussed further 
in Section 3.2. 

2.3. Line Fitting 

IRAF offers the widest range of possibilities: a 
line can be measured by summation of the pixel 
values (i.e., without any fitting procedure, Equa- 
tion 1), or by fitting a Gaussian function (Equa- 
tion 4), a Voigt profile, or a Lorentz profile. The 
task can measure both absorption and emission 
lines. Finally a deblending routine is also avail- 
able, where the user interactively marks the re- 
gion to investigate, the peaks of the blended lines, 
and so on. SPECTRE offers similar options, al- 
though it is more specifically designed for the clas- 
sical analysis of isolated absorption lines. ARES 
and EWDET use Gaussian profiles, performing 
fits similar to our Equation 4 and reporting the 
fit parameters in their output files. ARES, how- 
ever, appears to be able to fit a line with multiple 
Gaussian functions to improve the fit of each line 
profile, and the number of Gaussians actually used 
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is listed in the output file. 

Thus, the most widely used profile is a Gaus- 
sian function. Generally speaking, three parame- 
ters are free in the fit: the centroid Ao, the stan- 
dard deviation a or, equivalently, the full-width 
at half-maximum, and the depth of the line center 
(A in the notation of our Equation 4). Of these 
three, here we consider the FWHM. In an observed 
spectrum, the line profile is a convolution of the 
instrumental profile with an intrinsic line profile. 
At the resolution that is typically used for rou- 
tine abundance analysis (i.e., from R^ 20,000 to 
R/^ 50,000), the line profile is still dominated, for 
most lines, by the instrumental profile, which ei- 
ther is the same for all lines, or scales as the wave- 
length when multiple echelle orders are patched 
together. Therefore, in principle, there is no need 
to leave this parameter free for each individual line 
once the instrumental characteristics have been 
determined. 

Most existing packages, however, choose to 
leave this parameter free for each line, as is done 
in splot, SPECTRE, ARES, and EWDET. As an 
advantage, the derived FWHM can be used a pos- 
teriori to indicate whether a line has been prop- 
erly measured: lines with a FWHM significantly 
larger than average are most probably unresolved 
blends, while lines with a FWHM significantly 
smaller than average could be noise spikes, tel- 
luric features, or spectral defects. Another posi- 
tive side of this approach is that it compensates, 
to some degree, for saturation. When a line be- 
comes strong enough to just enter the saturation 
regime, it can still be mimicked with a slightly 
wider Gaussian, with a modest error in the EW 
measurement (of course, very strong lines cannot 
be correctly reproduced this way.) The main dis- 
advantage is allowing the FWHM to be completely 
free complicates the dcblcnding capabilities of an 
automated program when the noise level in the 
spectrum is appreciable. Most obviously, an ap- 
parently broad feature could be modeled by a 
single profile with a large FWHM or by multiple 
narrower profiles. Moreover, even when the num- 
ber of components in a blend is considered known, 
numerical degeneracies between the widths and 
strengths of the various overlapping features in- 
crease the uncertainties in the fitting parameters 
that are produced by the analysis. 



2.4. Uncertainty Evaluation 

Uncertainties in the measurement of EWs are 
often neglected when deriving chemical abun- 
dances. Often, having measurements of n lines 
for some chemical element, it is convenient to 
compute the mean and variance of the single-line 
abundances and to use oj\Jn — 1 as the uncer- 
tainty of the net abundance adopted for that par- 
ticular element. Using a good estimate of the EW 
uncertainty of each line to perform a weighted av- 
erage would give somewhat better results, but the 
dominant source of uncertainty in the adopted 
chemical abundances is often the uncertainty in 
the stellar atmospheric parameters. This explains 
at least in part why little effort is sometimes put 
into evaluating EW uncertainties (see also Sec- 
tion 2). 

Still, when writing a package that measures 
EWs, individual uncertainties can easily be esti- 
mated, and indeed they are of fundamental im- 
portance in particular cases, e.g., temporal vari- 
ations of the EW. In the case of splot, each of 
the pixels used in the fit is assigned its own stan- 
dard error through the construction of a statistical 
noise model based on Poisson photon noise (read- 
out noise is usually negligible for spectra with high 
S/N ratios). The error estimate, however, is op- 
tional and can only work if the user inputs an ap- 
propriate gain value if the spectral intensities arc 
recorded in some instrumental data-number units 
rather than in raw detected photoelectrons. Fur- 
thermore, if the input spectrum has been prenor- 
malized with a flat-field spectrum or to a provi- 
sional continuum, the numerical conversion from 
Poisson errors to intensity-unit errors will be a 
function of wavelength. In SPECTRE and ARES, 
no error computation routine is provided, for the 
reasons outlined above. EWDET, on the other 
hand, performs a full error computation that is de- 
scribed in detail by Ramirez et al. (2001). In short, 
each pixel used in the fit is assigned its own un- 
certainty estimate based on propagation of errors 
through the Gaussian fitting procedure 11 , includ- 
ing the continuum placement uncertainty, which is 
almost never considered in other packages. Pixel 



llr The Gaussian fit used by Ramirez et al. (2001) is different 
from our Equation 4, in that the formula by Ramirez et al. 
(2001) is strictly valid only when Ic=lj while our formula 
is valid in the general case, even if Iq ^1. 
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uncertainties are then summed in quadrature after 
multiplication by the pixel size. 

3. DAOSPEC 

DAOSPEC 12 is a code written in standard For- 
tran 77. The version that we have made generally 
available makes use of the cfitsio libraries 13 (Pence 
1999) for dealing with standard FITS-format spec- 
tra. IRAF format (.imh, .pix) files can also be 
used, and we have included an optional graphi- 
cal display 14 that allows the user to monitor the 
progress of the analysis. However, none of these 
should be considered an integral or essential part 
of the algorithms we discuss here. Practical and 
detailed help on how to obtain, install, and use 
DAOSPEC can be found in the DAOSPEC Cook- 
book 15 (Pancino & Stetson 2008), or by directly 
contacting the authors. 

We briefly describe here the general code struc- 
ture depicted as a flow chart in Figure 3. More de- 
tails on the most important steps are given in the 
following subsections. The reduction path consists 
of three main stages: (i) the input and prepara- 
tion phase, shown on the left side of Figure 3; (ii) 
the main iteration loop, center of Figure 3 and 
(Hi) the evaluation and output phase, right side of 
Figure 3. 

Three inputs, shown on the left side of Fig- 
ure 3, are needed: (i) a small set of configura- 
tion parameters, specifying among other things 
the wavelength limits of the spectral region of in- 
terest, a first guess at the FWHM, the order of 
the polynomial to be used for the continuum nor- 
malization, and other basic information required 
by the code; (ii) the reference line list, which con- 
tains the laboratory wavelengths of as many clean 
and unblended lines as possible, and may con- 
tain any additional textual information — such as 
atomic parameters and notes — that the user finds 
convenient; and (Hi) a spectrum in standard FITS 
or IRAF format, which must be binned linearly 
with wavelength. Items (i) and (ii) are specified 
once at the very beginning of the reduction pro- 



http://cadcwww.dao.nrc . ca/stetson/daospec/ 
http://www.bo . astro . it/~pancino/projects/daospec . html 
http: //heasarc .nasa. gov/docs/sof tware/f itsio/f itsio .html 
"SupcrMongo:" http : //www . astro . princeton . edu/~rhl/sm/ 
http : //cadcwww . dao . nrc . ca/stetson/daospec/ daospec . ps 
http://www.bo. astro . it/~panclno/docs/daospec.pdf 



cess; any number of individual spectra (Hi) may 
be then analysed one after another, at the user's 
pleasure. A more detailed description of this pro- 
cedure is given in the Cookbook. 

With the above inputs, the code performs an 
initial continuum fit with the help of Legcndrc 
polynomials of the order specified by the user, and 
performs the line finding procedure (see also Sec- 
tion 3.1). Then, a preliminary correlation of de- 
tected lines with the input list of laboratory wave- 
lengths is performed, to derive an initial radial- 
velocity estimate. 

Now the main loop starts (center of Figure 3). 
The detected lines are provisionally subtracted 
from the spectrum and the remaining residual 
spectrum is used to refine the continuum normal- 
ization (see also Section 3.2) by means of robust 
non-linear least squares. With the newly normal- 
ized spectrum, the individual line centroids and 
strengths arc then refined, along with a value for 
the FWHM which is either constant for all lines, 
or considered to be a linear function of wavelength 
(see also Section 3.3). This whole procedure is re- 
peated five times. (The number "five" was chosen 
arbitrarily and — so far — it appears to be adequate. 
It can be changed trivially if future experience in- 
dicates a need.) 

After final execution of the refinement loop, the 
fitted wavelengths of the individual lines are com- 
pared to the table of laboratory wavelengths. A 
simple outlier-clipping algorithm is used to deter- 
mine the final radial-velocity estimate for the star 
and, at the same time, to specify which of the de- 
tected lines lie sufficiently close to tabulated lab- 
oratory wavelengths to be regarded as legitimate 
detections of particular atomic transitions. Equiv- 
alent widths are computed from the fitted line pa- 
rameters, and uncertainties in those EWs are de- 
termined from the residual noise remaining in the 
spectrum within the wavelength range of the line 
profile (see also Section 3.3). 

The code output is a text file containing in its 
header the estimated FWHM at the center of the 
spectrum; the estimated radial velocity of the star 
with its standard error; the number of lines used to 
detemine the velocity; and the root-mean-square 
value of the pixel-by-pixcl flux residuals remain- 
ing in the spectrum after subtraction of all the 
fitted lines, expressed as a percentage of the fitted 
continuum flux. Following that, the file contains 
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a list of all the lines found by DAOSPEC having 
strengths greater than a user-specified minimum 
EW. Both their observed and inferred rest wave- 
lengths are reported, these latter computed on the 
basis of the estimated radial velocity. If a particu- 
lar detected line has been identified with an entry 
in the input line list, the laboratory wavelength of 
the line along with any additional textual infor- 
mation from the input line list is also reproduced 
in the output. 

Optionally, (right side of Figure 3), the code 
can be instructed to produce two FITS- or IRAF- 
format spectra, one containing the final estimated 
polynomial continuum model used to normalize 
the spectrum, and another containing the spec- 
trum of the flux residuals, of the flux residuals, 
after the spectrum has been divided by the fitted 
continuum and all the lines found by DAOSPEC 
have been fitted and digitally subtracted. 



3.1. Line finding, line identification and 
radial velocity 

Line finding in DAOSPEC is performed with 
the help of a tuned second-derivative filter to iden- 
tify local minima. The user-specified initial guess 
of the FWHM is used to optimize the distinction 
between valid lines and noise or continuum fea- 
tures 16 . During the main iteration loop (see Fig- 
ure 3) the initial FWHM estimate is refined at the 
same time as the line centroids and strengths, by 
the method of robust non-linear least squares, and 
provisional line detections are abandoned if refined 
EW estimates indicate they are too weak to be of 
interest. Note that, as mentioned above, the user 
may specify whether the FWHM is to be regarded 
as constant over the entire wavelength range of the 
spectrum, or as proportional to wavelength. 

Preliminary line identification is performed by 



This means that DAOSPEC is effective in separating noise 
from real spectral features as long as the FWHM of a line 
is adequately sampled, i.e., by at least 2 pixels. 
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a simple cross-correlation of the detected lines 
with the input list of laboratory wavelengths: that 
wavelength shift that corresponds to the great- 
est number of matching lines is adopted as in- 
dicating the preliminary radial-velocity estimate 
for the star. (The user can help make this pro- 
cess reliable by supplying as complete and appro- 
priate a line list as possible for the anticipated 
spectral class of the star, and for the presumably 
known wavelength range, resolution, and S/N ra- 
tio of the spectra.) After completion of the main 
refinement loop, these cross-identifications are re- 
considered on a line-by-line basis. A particular 
cross-identification of a detected line with a labo- 
ratory transition is regarded as valid if the radial 
velocity implied by that provisional identification 
agrees with the adopted final radial velocity of the 
star within a user-specified tolerance. Therefore, 
once the laboratory line list and the significance 
and identification tolerances have been specified 
by the user, both the line finding and identification 
procedures are completely automated, and the ob- 
served radial velocity of the star is one of the pro- 
gram outputs, not inputs. This has an impact on 
performace, as discussed in Section 2.1. 

The final radial-velocity measurement is pro- 
vided together with its line- by- line variance, (a), 
and the number of lines (n) actually used in 
its determination. These can in principle be 
combined to provide a final error estimate as 
a /\/n — 1. However, when measuring spectra 
that have not been specifically intended for high- 
precision radial-velocity measurements, special 
care must be taken. For example, if there are 
thermal changes between daytime lamp observa- 
tions and night-time star observations, or if the 
star image can be off-center in a finitely wide spec- 
trograph slit, there may well be systematic veloc- 
ity errors that are not reflected in the line-by-line 
velocity variance. 

There is the possibility of correcting for these 
effects if at least one of the telluric absorption 
bands 17 of O2 or H 2 is included in the observed 
spectral range. A very accurate list of rest wave- 
lengths of atmospheric absorption lines can be 



We remind the reader that while night-sky emission lines 
can be used to identify and correct thermal or other changes 
in spectrograph alignment, they cannot be used for correct- 
ing errors due to mis-centering of the star in the slit since 
the sky illuminates the slit uniformly. 



(a) Overestimated EW 


(b) Ideally Perfect EW 




^^^^^^^ 


(c) Effective Continuum 


(d) Improved EW 







Fig. 4. Effect of different choices of the contin- 
uum (see text for details). Panel (a): A strong 
spectral line is sketched as a Gaussian; small lines 
and spectral defects (small Gaussians) pollute the 
line; if the EW is measured with respect to the true 
continuum (horizontal line), it will be overesti- 
mated because the flux absorbed in the weak lines 
will be added to the flux absorbed in the strong 
line. Panel (b): as before, but now the shaded 
area represents the true amount of flux absorbed 
in the strong line alone, exclusive of the flux re- 
moved from the spectrum by the weak lines. Panel 
(c): here the strong spectral line has been fitted 
and removed, and the effective continuum (solid 
line) is obtained by balancing to zero the residu- 
als (dashed area), which include noise (not repre- 
sented here for clarity) and all the small pollut- 
ing features (small Gaussians); the effective con- 
tinuum (dotted line) is therefore lower than the 
true continuum. Panel (d): the EW is now es- 
timated with respect to the effective continuum, 
a depressed continuum that includes, statistically, 
the effect of typical local pollution by weak lines; 
the EW estimate is thus improved. 

obtained from the GE1SA 18 database (Jacquinet- 
Husson et al. 1999, 2005). A separate measure- 
ment can be performed using such a telluric line 



Gestion et Etude des Informations Spcctroscopiques 
Atmospheriques, http : //ara. lmd.polytechnique . f r/ 
htdocs-public/products/GEISA/HTML-GEISA/ 
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list as input and forcing DAOSPEC to explore a 
limited radial velocity range close to km s -1 . 
Our experience on red giants in w Centauri 19 (Pan- 
cino et al. 2007) and in open clusters 20 (Pancino 
et al., in preparation) shows that the radial veloc- 
ity shifts measured in the telluric absorption lines 
can be up to 1-3 km s _1 and, what is more im- 
portant, the associated uncertainty can be of order 

0. 5 km s _1 . However, when these measurements 
are used as corrections to the stellar radial veloci- 
ties, formal errors of 0.1-0.3 km s _1 can be easily 
achieved. 

3.2. Continuum Normalization 

What is generally understood as the contin- 
uum of a spectrum is the level of flux that one 
would have if only continuum-opacity sources — 

1. e, no discrete absorption lines or bands — were 
present: we will call this the true continuum. Tra- 
ditionally we seek the true continuum to anchor 
our fitting function when measuring the EW of a 
particular line (Equations 1 and 4). This is appro- 
priate when the spectra are not crowded and the 
line in question is not contaminated by neighbor- 
ing lines. Otherwise, true continuum determina- 
tions are particularly difficult and in most cases it 
is also very difficult to state the associated uncer- 
tainty. Moreover, we contend that when measur- 
ing lines in crowded spectra, what we really want 
to determine is not the true continuum level at all. 

In fact, we are trying to estimate how much 
flux one particular atomic transition is subtract- 
ing from the net output flux of the stellar atmo- 
sphere. But we need to separate that from the 
flux subtracted (or added) by other contaminat- 
ing effects, such as: (i) other atomic or molecular 
lines that are too weak or too close to our line 
to be distinguished; (ii) the noise, which for an 
infinite number of pixels should give an average 
zero contribution, but for a small number of pix- 
els could give a non-zero contribution; (Hi) detec- 
tor defects; and (iv) artifacts from the imperfect 
removal of cosmic rays, sky emission lines and tel- 
luric absorption features. Therefore, even if we 



Using ESO GIRAFFE data with R~20000 and S/N~50- 
100, and approximately 150 clean and unblended lines 
around the 6300A region. 

Using Calar Alto FOCES data with R~30000 and S/N~50- 
100, and a few hundred clean and unblended lines. 
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0.4 - 



0.2 - 



5536 553B 5540 5542 5544 

Wavelength (A) 

Fig. 5. — Artificial spectra used for continuum 
placement experiments, the clean spectrum (thin 
line) that includes only lines with EW between 
20 and 100 mA, and the full spectrum (thick 
line), that contains also several hundreds of small 
lines with EW<20 mA, the vast majority having 
EW<1 mA. 

were able to find the true continuum, when evalu- 
ating the area enclosed below the true continuum, 
we would not be measuring the flux of the intended 
transition alone, but its sum with all these other 
effects (Figure 4, Panels (a) and (b)). 

Rather than trying to find where the true con- 
tinuum lies, DAOSPEC looks for the effective con- 
tinuum, a depressed continuum that takes into 
account, in a statistical sense, the unrecognized 
flux deficits and excesses due to all the sources 
described above 21 . Under the assumption that 
the unidentified and unrelated opacity sources un- 
derneath the spectral line in question are statisti- 
cally similar to those lying at nearby wavelengths, 
DAOSPEC finds the effective local continuum by 
robustly balancing to zero the residual spectrum 
obtained after all the identified lines have been 
subtracted (Figure 4, Panel (c)). The area of a 
line measured with respect to the effective con- 



This is probably the reason why the effective continuum 
determined by DAOSPEC is sometimes lower than then 
continuum estimated by other programs, especially in the 
case of crowded spectra. Note however, that programs that 
rely on an automatic and global continuum placement pro- 
cedure, like EWDET, tend to be always in good agreement 
with DAOSPEC, placing the continuum lower than local 
continuum estimates. 
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Fig. 6. — Difference between the measured and 
true (input) EW on the artificial spectra. The left 
panels show DAOSPEC measurements, based on 
the effective continuum, and the right panels IRAF 
splot measurements based on the true continuum. 
Top panels show measurements on the clean spec- 
trum and bottom panels show measurements on 
the full spectrum. In each Panel, the dotted lines 
mark a perfect agreement (AEW=0), while the 
solid lines mark the median AEW. 

tinuum (Figure 4, Panel (d)) should be a more 
valid estimate of the correct EW of the particular 
atomic transition we are interested in, once the 
uncertainties of the procedure are properly taken 
into account. 

3.2.1. Experiment with Artificial Spectra 

We tested the above assertions with a simple ex- 
periment employing artificial spectra, created us- 
ing a list of more than 1000 atomic lines belong- 
ing to 50 species, in the wavelength range 5500- 
5600A. The atomic data were obtained from the 
VALD 22 database (Kupka et al. 1999). The most 
recent version of the spectral synthesis code origi- 
nally described in Spite (1967) was employed, to- 
gether with OSMARCS model atmospheres (Plez 
et al. 1992, and following updates), to predict the 
EW of each line for a metal-rich giant that should 
represent the typical difficult case: T e gf=4500 K, 
log 5 =1.0, « t =1.0 km s" 1 , and [M/H]=0.0, with 
no a-enhancement. The artificial spectra were 
generated as simply as possible, simulating lines 



http: //www . astro .uu. se/~vald/ 



with EW<100 mA with saturated Gaussians (Sec- 
tion 3.3), and chopping out stronger lines (dis- 
cussed in Section 3.3). A flat true continuum 
level was set to 1.0. The simulated spectral qual- 
ity is typical of a modern echelle spectrograph, 
FWHM=0.15 A (R ~ 40,000), S/N=35, a scale 
of 0.03 A/pixel. Two different spectra were gener- 
ated (Figure 5): a clean spectrum, where all lines 
smaller than 20 mA were omitted, and a full spec- 
trum including all lines, even the hundreds of tiny 
lines measuring as little as a fraction of a mA. 

The same list of 86 lines was then measured 
on both the clean spectrum and the full spectrum: 
once with DAOSPEC relative to its estimate of the 
effective continuum, and again with the task splot 
in IRAF relative to the true continuum, which is 
exactly known and fixed to 1.0 in this ideal case 
(see Figure 5). Results are presented in Figure 6. 
As anticipated, DAOSPEC's effective continuum 
is lower than IRAF's — by 1% in the clean spec- 
trum and by 2.5% in the full spectrum. Both dif- 
ferences are smaller than the corresponding root- 
mean-square residuals estimated by DAOSPEC, 
which are 2.4% per pixel for the clean spectrum 
and 2.9% per pixel for the full spectrum. Since 
the simulated noise was the same in both spec- 
tra, the slightly higher variance measured in the 
full spectrum is presumably due to the hundreds 
of additional small, unresolved lines that it con- 
tained. 

Both methods of measurement gave similar re- 
sults on the clean spectrum (top panels of Fig- 
ure 6): lines measured with respect to the effective 
continuum by DAOSPEC were underestimated by 
AEW=-1.60±0.95 mA while lines measured with 
respect to the true continuum by IRAF/ splot were 
overestimated by AEW=1.45±0.95 mA. Both dif- 
ferences are probably negligible compared to the 
uncertainties involved (~ 1.5<r). 

When considering the full spectrum however, 
apart from an overall increase in the spread of 
the measurements, we find that the DAOSPEC 
measurements relative to the effective continuum 
remain closer to the input values and more con- 
sistent with the measurements on the clean spec- 
trum than the IRAF/ splot measurements with re- 
spect to the known true continuum. Specifically, 
we find net shifts of (AEW) = -2.3 ± 5 mA 
for DAOSPEC/effective continuum and (AEW) = 
+6.0 ± 5.7 mA for splot/true continuum. We as- 
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sert, therefore, that as crowding increases, mea- 
surements relative to the true continuum — even 
when it is known perfectly — become unreliable be- 
cause of the likelihood of including opacity contri- 
butions from unrecognized and unrelated spectral 
lines in the EW measurement for the target line 23 . 

This is especially important to consider when 
one is measuring EW interactively with packages 
such as IRAF, MIDAS or SPECTRE. The proce- 
dure of readjusting the continuum locally around 
the line can produce a systematic bias if the user 
tries to ignore weak spectral lines in the continuum 
samples by favoring the highest points in the spec- 
trum or by selective placement of the continuum- 
sample cndpoints; this fails to allow for the possi- 
bility that flux may be removed by similar unin- 
teresting but nevertheless real absorption features 
occupying the same wavelength range as the line of 
interest. The depression of the continuum by such 
weak lines is, in fact, very difficult to estimate in 
an impersonal and repeatable way by eye when 
using only small regions around each line. How- 
ever, we believe that we can reduce the systematic 
errors due to unrecognized absorption features by 
iteratively using our model profiles to subtract all 
recognized spectral lines and then running the es- 
timated continuum through the middle of the re- 
maining flux spectrum, taking this to represent 
the typical local flux per pixel in the presence of 
all unmodeled opacity sources. 

3.3. EW Gaussian Fit 

When measuring EWs, DAOSPEC applies 
to crowded stellar spectra a concept also used 
in the analysis of crowded stellar images (e.g., 
DAOPHOT and ALLSTAR, Stetson 1987), i.e., 
the idea of fitting small groups of features simul- 
taneously. This concept is used also in the inter- 
active deblending routine of the IRAF task splot, 
but to our knowledge has not been implemented in 
an automatic program up to now. Besides allow- 
ing each absorption feature to be responsible for 
its own fraction of the flux missing from individual 
pixels, another advantage of this approach is that 



23 In reality, since perfect knowledge of the true continuum 
position in real — as opposed to simulated — spectra is im- 
possible, a person who attempts to estimate the true con- 
tinuum level to anchor EW measurements risks a system- 
atic error that is probably somewhere between the two lim- 
iting cases shown here. 



the fitting function of each feature now does not 
need to be truncated close to the feature center, 
but can extend outwards by several pixels (e.g., up 
to 6cr for a Gaussian function, where the intensity 
of a line is down by a factor e~ 18 <~ 10~ 8 from its 
maximum value), even in those regions where it 
overlaps one or more nearby features. Deblending 
as in DAOSPEC can therefore be more effective 
than in those programs that fit lines individually. 

The fitting function employed in DAOSPEC is 
based on a saturated Gaussian 24 defined as 



h(X) 



g(A) 
1+.9(A) 



(7) 



where g(X) is the Gaussian function (Equation 2). 
This particular formulation was chosen because 
deviations of h(X) from the original Gaussian 
shape are only important for strong lines: ac- 
cordingly, as g(A) becomes small h(X) tends to 
<?(A). Futhcrmore, this formula ensures that the 
residual flux at line center (~ 1 — h(X); cf. Eq. 3) 
approaches zero asymptotically as the EW grows 
without limit; if a pure Gaussian function had 
been retained as the model spectral line, there 
is no mathematical principle that prohibits the 
line profile from crossing over into negative fluxes 
(A > 1 in Eq. 2) as the line strength increases. 
Clearly, our adopted numerical formulation is not 
justified from physical first principles, but it is 
found empirically that it provides a useful exten- 
sion of the general method to moderately satu- 
rated lines. It cannot, of course, adequately re- 
produce lines that are strongly saturated. The 
exact EW at which the method is no longer valid 
clearly depends on the spectral resolution and S/N 
ratio, since the better a spectral line is defined, 
the more statistically significant its deviation from 
the model profile will be 25 . 



24 In the future, when higher resolution spectrographs will 
become routinely available, we plan to update DAOSPEC 
so that the user will be able to choose among different fit- 
ting profiles, including the Voigt profile, that would take 
care of lines with extended wings. It should be noted that 
these extended wings will substantially increase the anal- 
ysis times which may, we hope, be ameliorated by faster 
computers. 

25 Of course, not all lines deviate perceptibly from the 
adopted profile at exactly the same EW, since the line pro- 
file, and in particular the extended line wings, depend on 
the physical properties of the particular transition that pro- 
duces the line. 
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Fig. 7.— Comparison between DAOSPEC mea- 
surements and EWs by Rutten & van dcr Zalm 
(1984) on the UVES Solar spectrum for three 
different resolutions: R^IOOOOO (top panel), 
R-50000 (middle panel) and R-20000 (bottom 
panel). The EW at which the Gaussian approxi- 
mation adopted by DAOSPEC is no longer valid 
appears clearly as that EW where AEW starts to 
deviate from zero. 



To test how well our Gaussian approximation 
holds for real spectra, we have used the UVES So- 
lar spectrum 26 , which has an intrinsic resolution 
of about R~10 5 ; we have also degraded it to reso- 
lutions of R-50000 and R-20000. We have mea- 
sured EW with DAOSPEC using the list of clean 
Solar lines by Rutten & van der Zalm (1984). The 
difference between measurements by DAOSPEC 
and by Rutten & van der Zalm (1984) is shown 
in Figure 7 for the three Solar spectra. As can be 
seen, for the two high resolution cases DAOSPEC 
measurements begin to deviate from those of Rut- 
ten & van der Zalm (1984) at around 80-90 mA 
(for R-10 5 ) and 100-120 mA (for R-50000), 



26 http : //www . eso . org/observing/df o/ quality/UVES/pipel 
ine/Solar_spectrum . html 



while the agreement remains within the noise up 
to ^200 mA in the lowest resolution case. 

Another feature of DAOSPEC that helps de- 
blending in crowded spectra is the fact that the 
FWHM is consistent for all lines, instead of being 
left as a free parameter in the fit of each individ- 
ual line, as was already discussed in Section 2.1. 
Either the FWHM is forced to be the same for 
all lines, or it can scale with wavelength, which is 
more appropriate for cchclle high-resolution spec- 
tra that have been rebinned in wavelength. A first 
guess at the FWHM is input by the user and it 
is subsequently refined during the robust least- 
squares fits that take place in the main loop of 
DAOSPEC (central part of Figure 3). 

An additional advantage of the way DAOSPEC 
works again relates to the measurement of strong 
features in crowded spectra. For interactive rou- 
tines trying to fit each absorption line separately, 
the safe region is generally limited to the line core, 
which is most often approximately Gaussian. This 
leads to an underestimate of the EW of strong lines 
that have an important flux contribution in their 
(non-Gaussian) wings. DAOSPEC fits lines by 
minimizing the residuals over the entire line profile 
after all fitted lines have been subtracted: the total 
area of the residuals is forced to be as close to zero 
as possible (an oversubtracted core forced to bal- 
ance undersubtracted wings), and the total EW is 
still approximately correct to first order even if the 
line profile deviates from the adopted shape. This 
results from the fact that the whole line, includ- 
ing its extended wings, is included in the fitting 
region. Of course, as already noted, very strong 
lines cannot be properly measured this way. 

3.4. Uncertainties 

DAOSPEC provides three pieces of information 
that can be used to characterize the uncertainty of 
each EW measurement, described in the Sections 
below. 

3.4-1- Formal Standard Error 

The first is the formal standard error, cr(EW), 
obtained during the least-squares determination of 
the EW of each line. This is a purely empirical 
~ 68% confidence interval on the derived value 
of the line's EW, derived from the standard de- 
viation of the local flux residuals remaining af- 
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tcr the fitted continuum and all detected lines 
have been removed from the observed spectrum. 
Like EW, this has units of wavelength and re- 
sults from a standard least-squares propagation 
of errors that weights the residual of each indi- 
vidual pixel according to that pixel's contribution 
to the determination of the line's EW. The rel- 
ative error er(EW)/EW can be used as an effec- 
tive means of distinguishing good and bad lines, 
as shown in the top panel of Figure 8. Weaker 
lines have increasingly large relative errors, as ex- 
pected due to the presence of noise and random 
blending. When the error rises above, say, 20% 
the lines could be considered unsatisfactory and 
could be omitted from an abundance analysis 27 . 
Another way of using the relative error given by 
DAOSPEC would be to derive the abundance of 
one particular element not as a straight average of 
the single-line abundances, but as a weighted aver- 
age, using cr(EW)d(abundance)/d(EW) to define 
the weights. 

3.4-2. Quality Parameter 

The second is called the quality parameter Q, 
and is also associated with each individual EW 
measurement. For each line, the root-mean-squarc 
value of the intensity residuals is calculated for the 
same range of pixels as were employed in the pro- 
file fit for that line. The ratio of this root-mean- 
square residual to the root-mean-square intensity 
residual for the spectrum as a whole is that line's Q 
value. Thus, it tells whether the particular region 
where that line sits is more (Q > 1) or less (Q < 1) 
problematic than the average for the spectrum as 
a whole. 

The Q value for a given line can be higher than 
average if the line more strongly deviates from the 
Gaussian form than the average line, because in 
this case the local residuals will have a higher vari- 
ance than the rest of the spectrum. This is illus- 
trated in the middle panel of Figure 8, where the Q 
values of lines measured in star 2129 of the Open 
Cluster Cr 110 (from Pancino et., in preparation) 
are plotted as a function of EW. As can be seen, 
most lines above 200 mA have Q > 1, that is, the 
residuals near their positions are higher than the 
average residuals for other parts of the spectrum. 
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27 As a rule of thumb, lines with EW smaller than three times 
<t(EW) should certainly be rejected. 
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Fig. 8. — Behaviour of the relative error 
AEW/EW as a function of EW (top panel) and 
of the quality parameter Q as a function of EW 
(middle panel) and A (bottom panel). The plot 
reports DAOSPEC measurements for star 2129 in 
the Open Cluster Cr 110 (from Pancino et al., in 
preparation) . 



Another reason why Q could be higher than 
average is that the spectral quality (crowding, 
S/N ratio, defects) is not homogeneous along the 
whole wavelength interval. In the bottom panel 
of Figure 8, Q is plotted function of wave- 
length. As can be seen, lines at the blue end of 
the spectrum — where the S/N ratio is lower than 
in the rest of the spectrum — tend to have higher 
Q values than the rest. Also, the region around 
7600 A tends to have bad Q values due to the tel- 
luric absorption band of O2. In some cases (i.e., 
strongly varying S/N ratios) these effects could 
be alleviated by adopting different Q selections 
for different spectral regions. In some other cases 
(anomalous crowding, defects or molecular bands), 
it could be safer to cut out the whole region in A 
where anomalous Q values are clustered. 

As a general rule, therefore, a lintel in Q can 
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be defined, depending on the particular spectrum, 
to reject less reliable measurements. In the case 
of Figure 8, for instance, some value 1 < Qij m < 2 
could be adopted depending on the specific goals of 
the abundance analysis. However, please note that 
if a too strict selection in Q is adopted, some sys- 
tematic effects could be introduced in the analysis 
that are difficult to foresee. Perhaps safer would 
be simply to increase the adopted cr(EW) values 
for lines with moderately large values of Q in the 
ensuing abundance determinations. 

3.4-3. Standard deviation of Residuals 

The third parameter, which characterizes the 
quality of the solution as a whole and is given in 
the header of the output file, it is the standard 
deviation of the flux residuals, which can be con- 
structed as the r.m.s. relative (percentage) scat- 
ter in the effective continuum on a pixel-by-pixel 
basis. This has units of flux, and is mainly due 
to photon noise in the spectrum, of course, but 
substantial contributions can also come from un- 
recognized weak lines and molecular bands, which 
are generally more important in metal-rich, cool, 
and/or low-gravity stars, or from imperfections of 
the instrument or in the preliminary analysis, such 
as in the merging of echelle orders, the subtraction 
of night-sky emission lines, or the removal of CCD 
fringing effects. Note that since this index repre- 
sents a pixel-by-pixel residual variance about the 
adopted mean continuum, positive and negative 
EW errors should occur with similar frequency, 
and the overall error in the abundance of any given 
atomic species, due to these effects, should decline 
with increasing numbers of lines measured, expe- 
cially if the lines span a large wavelength interval. 

To evaluate the likely impact on the EW results 
due to the overall, systematic uncertainty of the 
continuum placement, a few considerations can be 
useful. First of all, the absolute uncertainty of the 
continuum level, in flux units, must be at least as 
large as the standard deviation of the flux residuals 
divided by the square root of the number of pixels: 



<t(Ic) > 



(t(1 pixel) 



Since it is often more convenient to work in rel- 
ative units, with a(l pixel)/ic representing the 
fractional standard deviation of the flux residuals, 



<t(I c ) o-(l pixel)// c 



Ic 



gives a lower limit to the relative uncertainty of 
the continuum placement. 

However, the N pixels have been used to derive 
a model continuum consisting of a polynomial of 
n th order, i.e., to determine n+ 1 free parameters. 
To be ultraconservative, the pixels that are cov- 
ered by spectral lines could also be entirely omit- 
ted from the uncertainty computation. There- 
fore, if ni is the total number of lines found by 
DAOSPEC and n is the polynomial order, the ul- 
traconservative estimate of the uncertainty of the 
global continuum placement procedure is 



o-(Ic) < cr(l pixel) /Jc 
Ic 



N-(n r FWHM) 
n+1 



(8) 



where we have supposed that each of the n + 1 
parameters of the continuum model is indepen- 
dently determined from at most pixels. For 
instance, in a modern spectrum containing 40,000 
pixels (iV) and 1,000 spectral lines (n;), having a 
FWHM of 10 pixels and a continuum shape cor- 
responding to a 20 t/l -order polynomial, and pro- 
ducing a residual spectrum with a noise level of 
1% (per pixel, root-mean-square), the overall un- 
certainty in the continuum placement cannot be 
as good as 1%/V40000 



0.005%, but neither is 
it likely to be as bad as 1%/ ^EEmm „ 
0.03%. The global impact of this continuum un- 
certainty on the EW results can be estimated by 
comparing a relative standard error calculated in 
this way to the typical biases illustrated in Fig. 2 
above. Alternatively, the spectrum, normalized 
with the continuum produced by DAOSPEC, can 
be multiplied by the corresponding 1 — AIc/Ic 
and by 1 + Ale I Ic and the two spectra thus ob- 
tained can be fed to DAOSPEC while switching 
off the continuum-fitting procedure (ORDER=- 
1). The difference between the two sets of EW 
and the original DAOSPEC result will give a more 
precise evaluation of the systematic uncertainty in 
the EWs due to the continuum placement proce- 
dure. 
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10 20 30 

Polynomial order 

Fig. 9. — An example that shows how the final 
output FWHM varies with the Legendre polyno- 
mial order (top Panel) and how the residuals of the 
subtracted spectrum (bottom Panel) vary with the 
Legendre polynomial order. 

3.5. Dependency on Input Parameters 

In this section we perform a few experiments 
to show how the results of DAOSPEC depend on 
the most relevant input parameters. A detailed 
description of the parameters can be found in the 
DAOSPEC Cookbook (Pancino & Stetson 2008), 
here we describe only the ones that have the most 
impact on the EW results. 

3.5.1. Polynomial Order for Continuum Fit 

As already mentioned in Sections 2, 2.2 and 3.2, 
overall continuum placement has a global and sys- 
tematic impact on the EW measurements. The 
only user-specified DAOSPEC input parameter 
that directly regulates it is the order of the Legen- 
dre polynomial used. This has little impact on the 
average continuum level because there are no other 
rejection or tweak parameters that alter the fit 
(see Section 3 for a description on how the contin- 
uum is fit). This is illustrated in Figure 9, where 
a sample spectrum has been fitted with different- 
order polynomials and the effect on the output 
FWHM and residuals of the subtracted spectrum 
has been shown. In general the FWHM changes 
are erratic, while the average residuals decrease 
when the polynomial order increases. 

To choose the appropriate order for the contin- 
uum fit, a good rule of thumb is to compare the 
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Input FWHM (pix) 

Fig. 10. — The three panels show (from top to 
bottom) how the number of lines identified in the 
input line list varies with the input FWHM, how 
the total number of lines found varies with the 
input FWHM and how the final output FWHM 
varies with the input FWHM. The dotted lines 
mark the correct FWHM=11.3 pixels. 

scale length of the global variations in the spectral 
shape to the total spectrum length. If the vari- 
ations are frequent and happen down to a scale 
length of 1/m of the total length of the spectrum, 
then the appropriate order to choose is m. If too 
low a value of m is adopted, there will be regions 
of the spectrum where the continuum is overesti- 
mated and other regions where it will be underes- 
timated, leading to an increase in the scatter of the 
resulting EWs (see also Section 4.1.2). For chem- 
ical species whose abundances rely on only a few 
lines, the mistake can be significant. On the other 
hand, provided that the spectrum contains many, 
many more than m valid pixels (which should nor- 
mally be the case) the main risk in adopting too 
large a value of m would be some increase in the 
reduction time. In case of doubt, we recommend 
that the optional output spectrum containing the 
fitted continuum be produced and overplotted on 
the original spectrum with, e.g., IRAF to check 
whether the order has been chosen appropriately 
(see also Figure 17). 

3.5.2. First-Guess FWHM and Scaling 

DAOSPEC needs a first guess for the FWHM 
(in pixels) to be able to initially distinguish and 
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Fig. 11. — A spectrum with a FWHM~5 pixels is 
measured with DAOSPEC using a FWHM which 
is underestimated by 10%, 20% and 30%. 

ignore features that are either too narrow or too 
broad to be likely real astrophysical spectral lines 
(see Section 3). Figure 10 shows what happens 
when the input FWHM is wrongly estimated by 
a large factor: while the number of detected lines 
that also appear in the input line list does not 
change very much, there is a significant change 
in the total number of lines found by the code; 
presumably this is because the input list contains 
primarily "interesting" — i.e., relatively strong — 
features, and it is the discovery of weak, mostly 
anonymous lines that is sensitive to the assumed 
FWHM. Since the working line list is limited to 
those lines detected initially, non-detection of real 
lines caused by initially specifying too large a 
FWHM is not fixed during the program's main 
iteration loop. However, spurious detections oc- 
casioned by the initial assumption of too small a 
value of the FWHM may be rejected later on, as 
the FWHM and the EWs are refined. We find that 
the output FWHM does not change significantly if 
the input FWHM is wrong by one pixel (in this 
case, ~10%), but starts to be over/underestimated 
perceptibly when the input FWHM is wrongly es- 
timated by more than that. Comparing the input 
and the output FWHM after a first run of the code 
can help in deciding whether the reduction should 
be reperformed with a different input FWHM. 

The characteristic FWHM of a spectrum is of- 
ten considered of secondary importance: we are 
not usually interested in its value except as it re- 



lates to the quality of the EWs that are the prin- 
cipal goal of the analysis. However, DAOSPEC 
could be used to estimate the FWHM of a spec- 
trum, which is extremely useful for spectral syn- 
thesis. While we have seen (Section 2) that a 
wrong average level of the continuum produces a 
systematic offset of the EW that does not depend 
strongly on EW, we show here that a wrong fitted 
value of the FWHM produces a systematic trend 
in the EW that depends on the EW itself. Fig- 
ure 11 illustrates the difference in the EW on lines 
measured on the clean spectrum described in Sec- 
tion 3.2.1 and Figure 5, when the input FWHM is 
kept fixed to a wrong value, i.e., underestimated by 
10%, 20% and 30% 28 . Therefore, if the FWHM to 
which DAOSPEC converges differs by more than 
approximately 10% from the user input FWHM, 
the measurement should be repeated until the in- 
put and output FWHM are sufficiently similar. 

One other important point to keep in mind 
is that in an echelle spectrum, where A/5 A is 
constant and the pixels are rebinned linear in 
wavelength, the FWHM changes from the blue 
to the red side of the spectrum. DAOSPEC al- 
lows the use of a constant FWHM or of a FWHM 
that scales with wavelength. Choosing a constant 
FWHM on an echelle spectrum would produce an 
overestimate of the FWHM in the blue and an 
underestimate on the red, with only the central 
EWs measured properly. This is illustrated in Fig- 
ure 12, where the echelle spectrum of star 2129 in 
Cr 110 (Pancino et al., in preparation), ranging 
from 5500 to 6900 A, is measured without scal- 
ing the FWHM and compared with the correct 
measurements where FWHM is scaled. As can be 
seen, the longer the spectrum, the bigger the mis- 
take, which can be of a few mA. When the EW 
differences are plotted as a function of wavelength 
(top panel) , we see a clear trend, with an increase 
of the spread on the sides, while when they are 
plotted as a function of EW (bottom panel) the 
only noticeable effect is an enormous increase in 
the spread. 



Please note that in this case, unlike in Figure 10, both 
the input and output FWHM are wrong, i.e., we force 
DAOSPEC to use the wrong values by setting FIX=1 
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Fig. 12. — EW difference on measurements done 
with constant FWHM and with FWHM scaled 
with A on an echelle spectrum of Star 2129 in 
Cr 110. In the top panel AEW is plotted versus 
wavelength, showing how EWs are overestimated 
in the blue and underestimated in the red. In the 
bottom panel AEW is plotted versus EW, where 
the mistake appears as an extremely large scatter. 



3.5.3. Residual Core Flux 

In real stellar spectra, the center of a line is 
never completely black because the effective sur- 
face of a stellar atmosphere is not at a temperature 
of absolute zero. Even strongly saturated lines, 
such as the H Q line (at 6563 A) in the spectrum of 
a low mass red giant, never reach zero flux even in 
their cores. Telluric absorption lines (i.e., the O2 
and H2O bands) or interstellar absorption lines 
(e.g., the NaD doublet) can, on the other hand, 
be much darker in their centers. Depending on 
the resolution and pixel size of a particular spec- 
trum, the maximum depth of strong lines varies 
a bit. For example, the spectrum of Figure 13 
has an H a line that reaches 20% of the continuum 
value in the line core. This effect is modeled in the 
DAOSPEC analysis by a further simple modifica- 
tion of Eq. (6) above: the model line profile now 
becomes 
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Fig. 13. — EW difference on measurements done 
with the correct Residual Core Flux value of 20% 
and the wrong value of 0%. The largest differences 
appear of course for strong lines. 



/(A) = 



3(A) 



l + g(X)/(l-p) 



(9) 



The model of the spectrum in the region of the 
line then becomes 

= I C ■ (1 - /(A)) - p Ic 

as 5(A) — > 00. The quantity /?, known as the 
"residual core flux," is specified by the user and is 
typically in the range 5-25%. It ensures that the 
profiles of strong lines saturate to a value equal to 
some fraction of the local continuum flux Ic rather 
than to a value of zero. A side-effect of this satura- 
tion model is that as lines become stronger within 
the saturation regime, they also become broader 
even for a fixed value of the FWHM parameter. 
This improves to some extent the capability of the 
model line profiles to mimic the actual profiles of 
real spectral features. 

3.6. Dependency on Spectral Quality 

To test the effect of the spectral quality on the 
EW measured by DAOSPEC, we created different 
versions of the full spectrum, our most crowded ar- 
tificial spectrum described in Section 3.2.1, with 
different S/N ratios, resolution (R=X/SX), and 
pixel sampling. 
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Fig. 14. — Difference on the EW measured with 
DAOSPEC and the "true" input EW of five artifi- 
cial spectra with R= 35000 and S/N ratios ranging 
from S/N= 300 (top panel) to S/N=10 (bottom 
panel). Continuous lines mark perfect agreement 
and dotted lines the a clipped average, also indi- 
cated on the top of each panel. 



Fig. 15. — Difference on the EW measured with 
DAOSPEC and the "true" input EW of five arti- 
ficial spectra with S/N=100 and resolution rang- 
ing from R=60000 (top panel) to R=5000 (bottom 
panel). Continuous lines mark perfect agreement 
and dotted lines the cr-clipped average, also indi- 
cated on the top of each panel. 



3.6.1. S/N Ratio 

The first experiment is shown in Figure 14, 
where the difference between the EW measured 
by DAOSPEC and the known input EW is plot- 
ted versus EW for five artificial R ~ 35000 spec- 
tra at various S/N ratios and with a sampling 
of approximately 5 pixels for each resolution cl- 
ement. Mean (cr-clipped) values with their line- 
by-line variance (a) are reported in Figure 14. As 
can be seen, there is no significant variation of the 
average AEW, which always remains close to zero, 
within the uncertainties. The variance a increases 
by approximately 2 mA from the S/N=300 to the 
S/N=10 case, as expected. Therefore, DAOSPEC 
gives good results even with a S/N ratio as low as 
10. This is due in part to the fact that the FWHM 



used by DAOSPEC is the same 29 for all lines: if 
the spectra are properly sampled 30 , DAOSPEC 
can usually tell a line from a noise feature. This 
result is somewhat in contradiction with a similar 
experiment performed by Sousa et al. (2006), and 
described in detail in Section 4.1.2. 

3.6.2. Resolution 

Resolution appears to be more critical than the 
S/N ratio as can be seen in Figure 15, where five 
spectra of fixed S/N=100 and similar sampling 
of approximately 5 pixels have been created at 
R ranging from 5000 to 60000. Again, the aver- 
age cr-clipped difference between DAOSPEC mea- 
surements and the "true" input EWs are always 



Scaled with wavelength in the case of echelle spectra, 
at least 2 pixels per resolution element. 
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consistent with zero within the spread. However, 
for R=5000, the total number of identified lines is 
very low. This, of course, is inherent to the in- 
trinsic quality of the spectra and not necessarily 
specific to the algorithms employed in DAOSPEC. 
Nevertheless, R=10000 can perhaps be considered 
a lower limit of validity of the EW method in gen- 
eral, and of the use of DAOSPEC in particular. 
We expect that stars with very high rotation ve- 
locities will behave as if they had a lower instru- 
mental resolution, as a consequence of the corre- 
spondingly broad line profiles. 

When resolution reaches R~40000, AEW gets 
definitely better not only in terms of the average, 
but also in terms of the spread, which at R=60000 
becomes ~3 mA. However, we must reiterate that 
our artificial spectra are simplistic, because all the 
simulated lines have pure Gaussian profiles. In 
real spectra, the effect of non-Gaussian wings on 
strong lines becomes important, as has been dis- 
cussed in Section 3.3 and shown in Figure 7. 

3.6.3. Pixel Sampling 

The last factor that may have an impact on 
EW measurements is the pixel sampling of the res- 
olution clement. The Nyquist-Shannon Theorem 
(Shannon 1949) states that the reconstruction of a 
continuous signal from its samples is possible if the 
signal is bandlimitcd and the sampling is greater 
than twice the signal bandwidth. In our case, the 
adequate reconstruction of a line profile is only 
possible if it is sampled by at least two pixels. 

Figure 16 shows the usual comparison between 
DAOSPEC measurements and "true" input EW 
on five artificial spectra with fixed R=35000 and 
S/N=100, and different sampling ranging from 
10 pixels (FWHM) to 1 pixel. As can be seen, 
all cases with FWHM > 2 pixels are comparable 
to each other, with the average AEW practically 
equal to zero and a ~ 8-9 mA. In the case of 
FWHM=1, the spread in AEW is naturally much 
higher, while the average AEW remains reason- 
ably close to zero. But a trend with EW seems 
present, in the sense that smaller lines tend to 
be slightly overestimated in a regime of under- 
sampling 31 . 

Since noise introduces variations at the pixel 
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31 This could be due to a larger uncertainty of the FWHM 
estimate in a regime of undersampling. 
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Fig. 16. — Difference on the EW measured with 
DAOSPEC and the "true" input EW of five arti- 
ficial spectra with S/N=100, R=35000 and pixel 
sampling ranging from FWHM=10 pix (top panel) 
to FWHM=1 pix (bottom panel). Continuous 
lines mark perfect agreement and dotted lines the 
a clipped average, also indicated on the top of each 
panel. 



level, whenever the sampling is close to the 
under-sampling limit (2 pixels), it is possible that 
DAOSPEC starts confusing noise with features, 
thus increasing the variance of measurements as 
in Figure 16 (bottom panel). Otherwise, the use 
of a correct first guess of the FWHM improves the 
separation of noise from real features. 

To summarize, DAOSPEC produces reliable re- 
sults for spectra that are adequately sampled, but 
the possibility of using DAOSPEC to measure EW 
also in spectra that are undcrsamplcd appears 
promising, although with a higher uncertainty (es- 
pecially for weak lines). 

3.7. Performance considerations 

As an automatic program, DAOSPEC is fast 
if compared to interactive or semi-automatic pro- 
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grams, especially when used in batch mode and 
without the graphical interface. However, given 
the large computational power of modern ma- 
chines, and given that it can be run 24 hours a 
day, 7 days a week without full-time supervision, 
DAOSPEC has been designed to give the most ac- 
curate results possible rather than the fastest pos- 
sible execution. In particular, the five big loops de- 
scribed in Figure 3 can take a long time in the in- 
terest of obtaining more reliable continuum place- 
ment, FWHM estimate, and line parameters. 

A few tests on different machines show that one 
spectrum can be measured in a short time, rang- 
ing from a few seconds to a few minutes, depend- 
ing on the spectral quality and other factors that 
are described in detailed in the following Sections. 
As an example, DAOSPEC took a few minutes 
to measure 2429 line on the FEROS spectrum of 
Section 4.1.2 (190 000 pixels long), a Solar-type 
star with a crowded spectrum and a complicated 
continuum shape. The above test was done on 
a Mac Pro workstation, with a Dual-Core Intel 
Xcon Processor (2x2.66 GHz - 4 cores) and 2GB of 
RAM memory, and the Mac OS X version 10.4.11 
(Tiger). 

3.7.1. Major factors 

The three major factors that have a large im- 
pact on the execution speed of DAOSPEC are the 
following: 

• The MXSPEC parameter, within the main 
code file (daospec.f ), specifies the maxi- 
mum number of pixels in the spectra that 
will be analyzed. The higher this value, the 
larger the amount of computer memory re- 
served. Setting MXSPEC to unnecessarily 
high values can cause the code to go unbear- 
ably slow or even to crash on some machines; 

• Setting the input FWHM to a severely 
wrong value can make the code go slower 
by a factor of 10 or even 100, as many exe- 
cutions of the inner iteration loops attempt 
to cope with the anomaly. We recommend 
deriving a robust estimate of FWHM before 
starting the computation, as suggested in 
the DAOSPEC Cookbook (Pancino & Stet- 
son 2008); 



• In general, the longer the spectrum (in pix- 
els) and the larger the number of lines in the 
spectrum (not in the laboratory line list), the 
longer it takes for DAOSPEC to converge. 
In complicated spectra, setting a larger value 
for the smallest interesting EW can reduce 
the number of loops and make the code go 
faster, although MXSPEC and FWHM have 
a much stronger impact. 

3.7.2. Secondary factors 

Other factors to take into consideration, al- 
though they have a smaller impact (see the DAO- 
SPEC cookbook, Pancino & Stetson 2008, for 
more details), are the following: 

• Adroit use of the Fortran compiler can opti- 
mize the code for faster execution; 

• Spectra with a simpler shapes, i.e., flatter 
continua, are faster to fit. Similarly, using a 
low continuum order (if appropriate) makes 
the computation faster, but it can also make 
execution significantly slower if used inap- 
propriately (see Section 4.1.2). Also, setting 
the continuum ORDER=-l (i.e., no contin- 
uum fit) can also speed up the computation 
if the user is satisfied that the spectrum has 
already been adequately rectified; 

• Using a FIXED FWHM makes the compu- 
tation a bit faster, but please see Figure 11 
for the risks associated with a wrong FWHM 
estimate; 

• Some very small improvement can also be 
obtained if the fixed FWHM option is used 
even for echelle spectra if they cover a very 
short wavelength interval, but at the expense 
of the EW accuracy (Section 3.5.2); 

• Using a smaller interval for the radial veloc- 
ity cross-correlation also makes the compu- 
tation faster. Leaving the search interval to 
its default value (±500 km s _1 ) is only ad- 
visable if there is no prior constraint on the 
approximate radial velocity of the objects; 

4. Comparisons 

In this Section we compare results produced 
by DAOSPEC with data and programs from the 
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current literature. We summarize and discuss all 
the papers that, to our knowledge, have made 
use of DAOSPEC in Section 4.1, including our 
own test with data from Pancino ct al. (2002). 
We also compare DAOSPEC with EWDET (Sec- 
tion 4.2) and with ARES (Section 4.3) and we fi- 
nally perform an abundance analysis of the Sun 
(Section 4.4). 

4.1. Literature tests on DAOSPEC 

DAOSPEC has been available to the astronom- 
ical community since 2002, when the first test ver- 
sions were circulated. Since then, it has evolved 
into the form presented here and in the Cookbook 
(Pancino & Stetson 2008), and it has been used 
and tested by some colleagues: a few authors used 
the code without mentioning any specific tests 
(Melendez et al. 2003; Pasquini et al. 2004; Dall 
ct al. 2005a,b, 2006; Pompeia et al. 2005; Zoccali 
et al. 2006; Lecureur et al. 2007). 

Other papers compare DAOSPEC measure- 
ments to manual measurements with IRAF, MI- 
DAS or other methods and find a good agreement, 
but they do not explicitly show the comparison 
(Pancino 2004; da Silva et al. 2005, 2006; Bar- 
buy et al. 2006, 2007; Letarte et al. 2006; Letarte 
2007; Alves-Brito et al. 2006). One paper used 
DAOSPEC to measure radial velocities and men- 
tions that a comparison with the results of fxcor 
within IRAF gives agreement within the uncer- 
tainties (Monaco et al. 2005). Another couple 
of papers publish the comparison of EWs mea- 
sured by DAOSPEC with manual measurements 
(Alves-Brito et al. 2005; Venn & Hill 2005); these 
are discussed in Section 4.1.1. 

To our knowledge, only two papers test DAO- 
SPEC extensively: Sousa et al. (2006) and Sousa 
ct al. (2007), which are discussed in detail in the 
following sections. 

4-1.1. Basic Comparisons 

Alves-Brito et al. (2005) used DAOSPEC to 
measure EW of Fe I and Fe II lines in five red 
giants in 47 Tuc. A comparison of IRAF inter- 
active measurements with DAOSPEC results on 
star 25 showed relatively good agreement, in the 
sense that the average AEW was smaller than the 
spread (cr=4.82 mA, see their Figure 3). But a 
slight trend with EW appeared, in the sense that 



for larger EW the AEW became larger. Accord- 
ing to the discussion in Section 3.5.2 above, an 
effect like this could, for instance, be due to a 
slightly (~10%) inappropriate input FWHM, but 
we do not know whether something like that is 
in operation here. In the end, the authors chose 
to adopt the IRAF measurements to derive their 
atmospheric parameters and iron abundances. 

Another, more favorable comparison was shown 
by Venn & Hill (2005), who plotted IRAF EW 
measurements by Shetrone et al. (2003) versus 
DAOSPEC, on GIRAFFE spectra (R-20000) of 
two stars in the Sculptor dwarf galaxy. They 
found good agreement (within 10%), with no sign 
of departures from the 1:1 relation for strong lines 
up to 200 mA. This is expected if one considers 
the example of Figure 7, where we show that the 
Gaussian approximation is more and more reli- 
able, even for strong lines, as resolution goes down 
from R~ 10 5 to R-20000. 

4-1.2. Detailed Comparisons 

Only two papers have performed detailed tests 
on DAOSPEC, namely Sousa et al. (2006) and 
Sousa et al. (2007). 

Sousa et al. (2006) used a synthetic (noise- 
less) model Solar spectrum of very high quality 
(R-120000 and S/N-300) to compare DAOSPEC 
and IRAF EW measurements. They found essen- 
tially perfect agreement in a red window (6000- 
6300A) with AEW=0.8±1.1 mA, based on 34 
lines, and fair agreement in a blue window (4400- 
4650A) with AEW=4.0±4.9 mA, based on 25 lines 
(see also Figure 20). This must of course be due 
to the higher crowding level and lower S/N ratio 
of the blue part of the spectrum. 

The synthetic spectrum was then degraded 
both in S/N ratio and resolution, and the DAOSPEC 
measurements were compared with each other. 
DAOSPEC appears to give very different average 
EWs and variance, by as much as AEW=15±20 mA, 
for the lowest resolution case (S/N~10 and 
R~ 12000). While some increase in the variance 
can be easily understood when varying S/N ra- 
tio or resolution, as can be seen in our own tests 
(Section 3.6), EW discrepancies and variance as 
large as those reported by Sousa et al. (2006) are 
difficult to understand, and indeed we do not find 
such behavior in our own tests (Sections 3.6.1 and 
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Fig. 17. — A typical FEROS spectrum of a Solar- 
type star (kindly provided by S. Sousa). Two 
Legendre DAOSPEC polynomials are overplotted, 
with an arbitrary vertical offset for clarity. As 
can be seen, the Legendre polynomial of 8 th order 
(lower continuous line) does not adequately repre- 
sent the spectral shape, while the 30 th order poly- 
nomial (upper continuous line) fits the spectrum 
better. 

3.6.2). 

When measuring their FEROS spectra, Sousa 
et al. (2006) encountered some problems. In 
particular, they found an enormous, unaccept- 
able spread in the resulting EWs, and they man- 
aged to obtain reasonable EWs only by cutting 
the spectrum into 100 A segments and running 
DAOSPEC manually on each small piece. The 
AEW from IRAF measurements then went down 
from 12.1±17.1 mA to 3.0± 4.7 mA (S. Sousa, 
private communication) but, of course, at the ex- 
pense of execution time and humanpower (several 
hours). They kindly provided us with some of 
their FEROS spectra (Figure 17) and we have re- 
peated their measurements. We found that using 
a different order for the continuum fit (30 instead 
of 8, Figure 17) and a different FWHM (14 instead 
of 5, Section 3.5.2, Figures 10 and 11) gave much 
better EWs and decreased the execution time by 
a factor of 50, roughly. We also tried cutting the 
spectrum into short pieces, both as a consistency 
check and to test the execution times, but we used 
shell scripts (Cookbook, Pancino & Stetson 2008) 
to run DAOSPEC automatically — in 10 minutes, 
total time (see Section 3.7) — on the various pieces: 
we obtained AEW=-4.1±4.3 mA when using the 



full spectrum, and AEW=-6.5±4.4 mA when the 
spectrum was cut into 100 A pieces. 

The paper that introduced ARES (Sousa et al. 
2007) was the second to perform a detailed check 
on DAOSPEC, using the same datasets as Sousa 
et al. (2006) and the same DAOSPEC configura- 
tion parameters. ARES is based on the IRAF 
task splot, and therefore the first comparisons 
made were between ARES and IRAF, and between 
DAOSPEC and IRAF. The results obtained with 
ARES were more similar to IRAF than the ones 
with DAOSPEC, supporting the conclusion that 
ARES is a very well designed extension of splot. 
We have seen, however, that the most important 
factor in these comparisons can be the way the 
continuum is chosen. In case of crowded spectra, 
we have claimed that the algorithm employed by 
DAOSPEC can give better results (Section 3.2.1), 
but since both IRAF and ARES are highly cus- 
tomizable in terms of continuum placement, we 
do not doubt experienced and careful users can 
obtain good results with those algorithms. 

To summarize, an appropriate choice of the con- 
figuration parameters is crucial to obtain good re- 
sults with DAOSPEC. The Cookbook (Pancino 
& Stetson 2008) provides practical and objective 
methods for finding the best values for these pa- 
rameters, as well as the discussions and tests pre- 
sented in Section 3.5 here. 

4-1.3. Red Giants in lo Centauri 

The dataset of EW measurements that Pancino 
et al. (2002) obtained with the IRAF task splot to 
derive abundances for six red giant stars in u Cen 
constitutes a good testbed for DAOSPEC. The 
full data description can be found in the original 
paper; in short, the six spectra were taken with 
UVES at the Very Large Telescope in Paranal, 
Chile, with R ~45000 and S/N~ 100-150 per res- 
olution element, covering the range 5250-6920A. 
Stellar metallicities range from [Fe/H]=-0.49 to 
-1.20, with temperatures around 4000 K and grav- 
ities of about 1 dex. The input line list contains 
230 features of various elements, although only 
[Fe/H], [Ca/Fc], [Si/Fe] and [Cu/Fe] were pub- 
lished by Pancino et al. (2002). 

We remeasured these spectra with DAOSPEC 
and compared the results (Figure 18). A to- 
tal of 1150 lines were used in the comparison. 
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Fig. 18. — Comparison of the original measure- 
ments from Pancino et al. (2002), obtained with 
IRAF (Y axis), and the measurements obtained 
here with DAOSPEC with the same line list and 
on the same spectra (X axis). Perfect agreement 
is marked with a dotted line. 

We found a very good average agreement, with 
DAOSPEC measurements marginally smaller, by 
AEW=-1.3±10.3 mA. When considering the six 
stars separately, we found differences ranging 
from AEW=-3.7±10.7 mA, for star WFI 222068, 
which is the most metal rich of the sample, to 
AEW=1.1±7.1 mA, for star WFI 618854, which 
is the most metal-poor of the sample. No trend 
with EW is apparent. 

The agreement appears satisfactory within the 
uncertainties, expecially in the light of the tests 
performed in Section 3.2.1, where we show again 
that an agreement between DAOSPEC and IRAF 
measurements gets naturally worse as metallicity 
(and line crowding) increases. 

4.2. DAOSPEC vs. EWDET 

EWDET (Ramirez et al. 2001, see also Sec- 
tion 2) was obtained by courtesy of S. Ramirez. 
It came with a test spectrum of a moderately 
crowded red giant in M 71, covering the range 
7900-8000 A, with R~30000 and a S/N ratio of at 
least 100 everywhere. We used this spectrum with 
the default configuration file ewdet.inp provided 
with the code, to measure the EWs of 70 lines. 
Ten additional lines were found, but EWDET did 
not report an EW for any of them because the 
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Fig. 19.— Difference between DAOSPEC and 
EWDET EWs measured on a test spectrum of a 
giant in M 71 provided with EWDET. AEW is 
plotted versus EW (top panel) and versus wave- 
length (bottom panel). The average difference is 
AEW=-1.2±11.7 mA. 

Gaussian fit did not converge. All the lines found 
by EWDET were used as the input "laboratory" 
line list for DAOSPEC, which we then used to 
obtain EWs from the same spectrum. It is per- 
haps worth stressing here that the input line list 
plays no part in the finding of candidate spectral 
lines by DAOSPEC; it is only after features have 
been detected that tentative identifications with 
features in the input list are sought. There is no 
attempt to "force" the detection of features in the 
spectrum at wavelengths specified by the input 
laboratory list. In the present case, DAOSPEC 
was able to (independently) find and measure all 
the lines that EWDET had found, including the 
10 that EWDET had subsequently discarded. No 
apparent defect was found on those 10 lines in a 
visual inspection of the spectrum. 

Figure 19 plots the difference between the 
DAOSPEC and EWDET EWs versus EW (top 
panel) and versus wavelength (bottom panel). 
The average difference is AEW=-T.2A with a 
variance of 11.7 mA. While the two sets of mea- 
surements appear in good agreement, the spread 
is slightly higher than expected, i.e., higher than 
that found in the comparison of Section 4.1.3 be- 
tween DAOSPEC and hand-made IRAF measure- 
ments. A trend of increasing spread with increas- 
ing EW might be present, while no obvious trend 
with wavelength is seen. 
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A comparison of the model continua adopted 
by the two programs shows a an overall systematic 
difference of 1.3% (DAOSPEC continuum lower), 
with a variance of 0.6% around this mean offset. 
For comparison, the residual spectrum produced 
by DAOSPEC has a pixel-to-pixel flux variance of 
2%. On the surface, this case appears to be similar 
to that discussed in Section 3.2 above. Such a 
discrepancy in the continuum levels could be the 
cause of the small AEW offset found between the 
two codes (see Section 2, Figure 2). 

On average, the standard errors estimated by 
DAOSPEC are larger by 0.7±1.5 mA than those 
reported by EWDET, even though the latter also 
includes the uncertainty due to the continuum 
placement and the former does not. In any case, 
given the large spread in AEW seen in Figure 19, 
both error estimates appear a bit small, indicat- 
ing that some other unidentified source of uncer- 
tainty might be present. If we estimate an error 
budget, including the average errors by EWDET 
(-3 mA) and DAOSPEC (-4 mA) and an er- 
ror due to the continuum placement as estimated 
roughly from Figure 2 (<~7 mA), we account for a 
spread of ^9 mA, i.e., the missing source of un- 
certainty must be of about ~8 mA 32 . This might 
suggest that the EWDET continuum placement 
uncertainty (Ramirez et al. 2001, and Section 2.2 
here) might be underestimated. 

Finally the average difference between the 
FWHM found by EWDET for each line, and the 
FWHM found by DAOSPEC (scaled with wave- 
length) is AFWHM=0.001±0.076 A, and the av- 
erage radial velocity difference between the two 
sets, in the sense EWDET minus DAOSPEC, of 
measurements is very small, Aiv=0.1±0.6 km s _1 . 

Summarizing, the comparison can be consid- 
ered satisfactory once all the sources of uncer- 
tainty are properly taken into account. The only 
minor disadvantage of EWDET is the fact that it 
has been written for personal use and it requires 
knowledge of Fortran to manually adapt some rou- 
tines to meet the needs of each set of spectra, in- 
cluding naming conventions and so on. 
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Fig. 20.— Difference between DAOSPEC and 
ARES EWs measured on a test spectrum of the 
Sun provided with ARES. AEW is plotted ver- 
sus EW (top panel) and versus wavelength (bot- 
tom panel). The average difference is AEW=- 
1.1±3.7 mA. 



4.3. DAOSPEC vs. ARES 

ARES (Sousa et al. 2007) has been obtained 
from the ARES webpage 33 with, inter alia, a test 
spectrum of the Sun obtained with HARPS from 
observations of Ganymede 34 . Similarly to what 
we have done with EWDET, wc ran ARES on the 
test spectrum and used its output as an input line 
list for DAOSPEC. 

The average differences of key parameters, in 
the sense of DAOSPEC minus ARES, can be 
summarized as follows: AEW=-1.1±3.7 mA, 
AFWHM=0.01±0.05A 35 , Aw r =-0.002±0.126km s"\ 
based on 98 lines in common. No error estimate is 
provided by ARES. At first glance, all these values 
appear in very good agreement within the uncer- 
tainties, even better than the comparison made 
with EWDET in Section 4.2. This is especially 



32 One thing to note here is that the spectral range chosen 
includes a relatively strong H2O telluric absorption band, 
which could cause some of the observed scatter. 

33 http:/ /www. astro. up.pt/~sousasag/ares/ 

34 http: / /www.ls. eso.org/lasilla/sciops / 3p6/harps / monitorin 
g/sun.html 

35 This difference is 100 times larger than in the case of the 
EWDET comparison. However, it corresponds to a ~8% 
relative difference, which, according to our Figure 11, im- 
plies a AEW of a few mA at most, compatible with the 
quoted AEW=-1.1±3.7 mA. 
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s Tablc 1 
Abundance Analysis of the Sun 
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true when considering the spread in AEW, which 
is 11.7 mA in the comparison with EWDET and 
only 3.7 mA in the comparison with ARES. The 
very good agreement must of course be largely due 
to the fine quality of the test spectrum, which has 
R~45000 and S/N~350. Figure 20 confirms good 
agreement with no trends with wavelength or EW 
in the differences, except for a possible problem 
in the bluest and most crowded part of the spec- 
trum. A last comparison was made on the number 
of lines found. The authors do not mention how 
many lines were found and/or identified by each 
code, but state that ARES finds more lines than 
DAOSPEC. If we compare the Solar spectra taken 
from Ganymede, we find that ARES identifies 101 
lines, and DAOSPEC identifies 100. 

ARES and DAOSPEC represent two very dif- 
ferent ways of approaching the problem of measur- 
ing EWs. ARES closely follows IRAF, including a 
major IRAF feature, namely the possibility to cus- 
tomize the continuum fitting procedure. Because 
of the way the continuum is fit, ARES is faster 
than DAOSPEC, although maybe a bit longer to 
configure. ARES takes of the order of seconds for 
each spectrum, while DAOSPEC may take from 
a few seconds to a few minutes, depending on the 
spectrum characteristics. Finally ARES gives no 
error estimate or radial velocity; indeed, the radial 
velocity is one of the necessary inputs, not outputs 
of the code. 

Nevertheless, in spite of the different continuum 
placement philosophies, ARES and DAOSPEC 
give entirely comparable measurements, within 
the uncertainties. 

4.4. Abundance analysis of the Sun 

As a final test, DAOSPEC was used on the So- 
lar spectrum obtained with HARPS (Section 4.3) 



to derive iron abundances for the Sun. The re- 
sults have been compared, using the same mod- 
els and abundance calculation code, to the abun- 
dances obtained with the EWs measured by Moore 
ct al. (1966) and Rutten & van der Zalm (1984). 

To measure EWs with DAOSPEC, we created 
a line list containing all lines in common between 
Moore et al. (1966) and Rutten & van der Zalm 
(1984). This line list was fed to DAOSPEC and, 
for homogeneity, it was also used to derive the So- 
lar abundance with the original Moore et al. (1966) 
and Rutten & van der Zalm (1984) measurements. 

We used the atmospheric models by Edvards- 
son ct al. (1993) and the latest version of the 
abundance calculation code originally published 
by Spite (1967). For sake of homogeneity, the 
atomic parameters (including log <?/) were taken 
from the line list of Rutten & van der Zalm (1984). 
In this way, the only difference among the three 
analyses comes from the EWs. The Solar temper- 
ature was kept fixed at 5780 K; gravity was al- 
lowed to vary between log 3=4.4 and 4.5, to allow 
for micro-adjustments of the Fe I and Fe II ioniza- 
tion equilibrium; the micro-turbulent velocity was 
kept as a free parameter and the best value was 
chosen as the one that minimized the slope of the 
EW vs. [Fe/H] relation. 

The results of our analysis are shown in Table 1, 
where it can clearly be seen that only negligible 
variations in log g and very small variations of 
the microturbulence (v t ) were necessary to obtain 
abundances that are quite compatible with each 
other, and with the Solar values. 

5. Conclusions 

While we have argued that one of the advan- 
tages of DAOSPEC is that it is minimally inter- 
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active, and hence fast and reproducible, we must 
confess that there still remain some user-specified 
reduction options whose best values are not en- 
tirely obvious. Examples of these are the polyno- 
mial order to be used in the continuum fits, the 
initial estimate of the FWHM value appropriate 
to a given spectrum, the range of plausible radial- 
velocity values to consider, and the minimum EW 
that is considered to be of interest. In some sense, 
however, we feel that this evident weakness is also 
a not-so-evident strength: the program runs suffi- 
ciently fast that the user is easily able to alter each 
of the various input parameters over some reason- 
able range of values and see directly how the astro- 
physical results are changed when all other inputs 
are held strictly fixed. We have also explored in 
detail the behavior of DAOSPEC in several cases 
of interest, performing a variety of tests that show 
the range of applicability of DAOSPEC in terms of 
data quality, and should also help in the search for 
the best set of configuration parameters for each 
case. Novice (and experienced!) users are posi- 
tively encouraged to carry out such experiments. 
Similarly, feedback from the users is of immense 
help, allowing us to refine the code in the direc- 
tions of major interest for the scientific commu- 
nity. 
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